Universality of liquid dynamics 
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We investigate the origin of the Stokes-Einstein relation in liquids. The hard-sphere dynamics is 
analyzed using a new measure of structural relaxation - the minimum Euclidean distance between 
configurations of particles. It is shown that the universal relation between the structural relaxation 
and diffusion in liquids is caused by the existence of one dominating length scale imposed by the 
structural correlations and associated with de Gennes narrowing. We demonstrate that this relation 
can be described by a model of independent random walkers under the single-occupancy constraint. 
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Much effort in the studies of liquids has been fo- 
cused on understanding the relationship between the 
crgodicity-restoring structural relaxation and diffusive 
motions of individual particles [l|, |2j. It has been rec- 
ognized that, at least in simple liquids, the self-diffusion 
coefficient D can be derived from static correlation func- 
tions without any knowledge of the interparticle forces 
0, 0. On the other hand, kinetic theories 0, Q have 
been successful in expressing time-dependent correlation 
functions describing the relaxation dynamics in liquids 
in terms of structural correlations. These results sug- 
gest a purely structural mechanism linking the two as- 
pects of liquid dynamics. The existence of such link can 
also be inferred from the Stokes-Einstein |7j and Stokes- 
Einstein-Debye || relations which can be combined in 
the following relation [9|, [l0( : 
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where r is the characteristic time of structural relaxation, 
and a is the size of the diffusing particle. We note that in 
dense fluids, a can be interpreted as the close-neighbour 
distance which represents the dominating length scale of 
the structural correlations. 

The current interest in the mechanism of the Stokes- 
Einstein relation is motivated by the recent finding that 
this relation breaks in strongly supercooled fragile glass- 
forming liquids 0, flil HH ] . A number of models have 
been proposed to interpret this observation 0, H3 |. 
However, we must first understand why these relations, 
derived from purely hydrodynamic postulates, demon- 
strate such remarkably universal success in describing the 
atomic-scale dynamics in the normal liquid state. 

In this Letter, we analyze the origin of the Stokes- 
Einstein dynamics in liquids. The dynamics of hard 
spheres (HS) is examined using a new measure of struc- 
tural relaxation based on the assessment of minimum dis- 
tance between the system's configurations. We present 
evidence that the invariance of the relaxation rate rela- 
tive to the diffusion rate in the liquid state is linked to 



the existence of a single structure-induced length scale 
dominating the liquid dynamics which is manifested by 
de Gennes narrowing |15| , and it also implies Gaussianity 
of the diffusion process. We show that this dynamics can 
be described by a simple model of independent random 
walkers confined to single-occupancy cells ^| 0] . 

The dynamics of a system of identical particles is com- 
monly described in terms of the van Hove density corre- 
lation function G(r,t) [lj. Its self-part, G s (r,t) describes 
the evolution of the density distribution for a tagged par- 
ticle. However, because the particles are indistinguish- 
able, the evolution of a tagged particle becomes irrele- 
vant for the structural relaxation as soon as the particle 
traverses the nearest neighbour distance. 

For a system of N identical particles, there exists N\ 
physically indistinguishable permutations of each config- 
uration. Let r N (t) = {r\;(i)} be the system's configura- 
tion at a time moment t. We define a correlation func- 
tion: 



1 / N 

G n {r,t) = -(Y j 5[r-\vm-Mi)( t )\ 



(2) 



\y( i)\ = P{i}, where P is the permutation procedure 
[18| minimizing the second moment: 
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This quantity, to be referred to as dissimilarity, repre- 
sents the minimum square Euclidean distance between 
the two configurations. In general, this minimizing re- 
labeling procedure (d-mapping) can be used to evaluate 
the dissimilarity between any two configurations: 
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We use this approach to analyze a molecular dynamics 
model of 27000 identical HS particles. Throughout this 
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study, the length is measured in terms of p 1 ^ 3 , p being 
the number density. 

For high packing fractions and short time intervals 
where the number of relabellings is negligible, G n (r,t) 
would behave as G s (r,t). Therefore, the former, like the 
latter, is expected to have a Gaussian form. However, 
unlike G s (r,t) the second moment of which - the mean 
square displacement (MSD) - grows linearly with time, 
G n (r, t) converges to a limit distribution with a finite sec- 
ond moment d(t) — > di as the two configurations become 
uncorrelated for t — > oo . The variation of di in the HS 
system is plotted in Fig. as a function of the packing 
fraction tp. 

For a selected reference configuration, the described 
permutation P uniquely maps the system's configura- 
tion space onto a compact region possessing the following 
properties: 

(i) For N — > oo almost all its points are infinitely close to 
the surface of a sphere of square radius Ndi 01 centered 
at the reference configuration (d-sphere). Therefore, any 
configuration uncorrelated with the reference configura- 
tion is expected to be separated from the latter by d — di . 

(ii) di separating two uncorrelated configurations is a 
static property entirely determined by the density fluc- 
tuations in these configurations. In the HS system, it 
increases monotonously with decreasing packing fraction 
if, see Fig. ^ The impact of the structural fluctuations 
on di can be assessed by using a regular lattice as a ref- 
erence configuration. It is shown in Fig. ^ that for the 
entire range of tp, di for the HS system calculated with 
respect to the fee lattice is reduced compared with di 
calculated for two uncorrelated HS fluid configurations. 

(iii) In the case of the perfect gas (PG), the volume of 
the d-sphere scales as the total number of distinguishable 
configurations. To test this, we consider a PG model 

;o si ngle -' 

mo. 

the fee configuration is equal to the inertia moment of 
the cell. The single-occupancy constraint reduces the 
total number of PG configurations by the factor e N [l7j : 
respectively, di is expected to be reduced by the factor 
e 2 / 3 . As shown in Fig.^ this is in good agreement with 
the simulation results. 

We now prove that configurations r N (0) and r N (t) are 
uncorrelated if d(t) — di. Consider d-mapping of r N (t) 
onto the d-sphere of a reference configuration c N such 
that all its particles have identical environments. Rela- 
tive to a pole of a multidimensional sphere, almost all the 
points of the sphere are infinitely close to the respective 
equator. Therefore, if d(t) = di, r N (i) would be mapped 
on the equator assuming that r JV (0) is mapped as a pole. 
The orthogonality condition follows: 



where particles are confined to single-occupancy Wigner- 
Seitz cells of the fee lattice O |T3 . Its di relative to 
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FIG. 1: Dots, di for two uncorrelated HS fluid configurations; 
triangles, di for a HS fluid configuration relative to the fee 
reference configuration. The arrow indicates di for perfect 
gas constrained to the single-occupancy fee Wigner-Seitz cells 
Hi El multiplied by e 2/3 . 



where p(i) and q(i) are permutations d-mapping r Af (0) 
and r N (t) onto c N . Covariance of two stochastic vari- 
ables is a measure of their linear correlation; therefore, 
(5) implies that r N (0) and r N (t) are uncorrelated (no- 
tice that d-mapping can induce correlation between these 
configurations if c N includes any density fluctuations). 
Fig. Elshows the evolution of the l.h.s. of Eq. (5) in the 
HS liquid for tp = 0.48 as a function of d(t)/di, the refer- 
ence configuration being the fee lattice. 

Thus, d(t) = di is a necessary and sufficient condition 
for the two configurations separated by d(t) to be un- 
correlated. d(t) can therefore be used as a convenient 
indicator of the system's approach to its ergodic equilib- 
rium. We propose to measure the structural relaxation 
in a system of particles in terms of the following dimen- 
sionless quantity: 



$ d (t) = [di - d(t)]/di 



(6) 
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Structural relaxation in liquids is commonly assessed in 
terms of the intermediate scattering function F(Q mi t), 
Q m being the position of the main maximum of the struc- 
ture factor S(Q m ) = F(Q m ,0) USUI. For liquid den- 
sities, Q m represents the length scale of the slowest re- 
laxing structural correlation (De Gennes narrowin g | 15| ) . 
This effect, however, disappears at low densities flLl2Lll5l| . 
An important advantage of as a measure of relax- 

ation is that it is not associated with a particular length 
scale. Also, unlike the energy-based measure of ergodic 
convergence [^j, its calculation doesn't require knowl- 
edge of the interparticle interaction. Another advantage 
of 3>d(t) is that, being based on a single-particle correla- 
tion function, it can provide insight into the link between 
the individual particle dynamics and the structural relax- 
ation. 
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FIG. 2: The evolution of the scalar product in the l.h.s. of 
Eq. (5), as a function of d(t)/di, calculated for the HS liquid 
for (p = 0.48. The reference configuration is the fee lattice. 
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FIG. 4: S(Q)/Q 2 for the HS fluid representing the Q- 
dependence of the relaxation time of the intermediate scat- 
tering function F(Q,t) |24l ]. Solid line, ip = 0.2; dashed line, 
p> — 0.3; chain-dashed line, <p = 0.35. 
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FIG. 3: Relaxation in the HS system as a function of MSD. 
Dots, filled boxes and filled stars: &d{t) for ip = 0.48, ip = 0.35 
and ip = 0.2, respectively. Open circles, open boxes and open 
stars: F(Q m ,t)/S(Q m ) for = 0.48, <p = 0.35 and p = 0.2, 
respectively. Crossed circles, <&d(t) for perfect gas. Dashed 
line, &d(t) calculated according to Eq. (7) 



In Fig.|3|we analyze the relation between the structural 
relaxation and diffusion in the HS system. The structural 
relaxation is measured in terms of &d(t) and presented 
as a function of MSD. For two representative values of 
the packing fraction, ip = 0.35 and ip — 0.48 which span 
the liquid domain, the evolution of 3>d(t) is compared 
with the respective evolution of F(Q m ,t)/S(Q m ). It is 
evident that for these two values of ip, the decay of $d(t) 
is exponential, and it agrees quite well with the respective 
behaviour of F(Q m ,t)/ S(Q m ). It is also evident that, 
within the liquid domain, the rate of structural relaxation 



relative to the diffusion rate is invariant with respect to 
the variation of ip. 

Fig. |3| shows that the universal relation between the 
structural relaxation and diffusion observed within the 
liquid domain breaks for lower densities. At p> = 0.2, 
we find a distinctly different regime. The structural re- 
laxation relative to the diffusion as measured by 
becomes retarded as compared with the universal pat- 
tern observed for the liquid. Notice that the decay of 
F(Q m ,t)/S(Q m ) for ip = 0.2 agrees quite well with the 
exponential liquid pattern. One can also see that $>d(t) 
exhibits a strong stretching as compared with the ex- 
ponential decay. Both the retardation and the degree 
of stretching progressively increase with decreasing <p. 
These results can be compared with an earlier simula- 
tion study of the HS viscosity [2l| where it was found 
that the Stokes-Einstein relation breaks around (p = 0.15. 
We note that the Enskog theory j^l predicts violation of 
the SE relation in the low-density gas limit, with Drj di- 
verging as (p^ 1 - 

These results thus pose us with the following questions: 
what is the origin of the universality of the liquid relax- 
ation dynamics and why does it break in the low-density 
domain? 

The observation that, within the liquid domain of ip, 
&d(t) decays exponentially relative to MSD implies that 
the liquid relaxation is controlled by a single length scale, 
and the close agreement between the decay of $>d(t) and 
the decay of the respective F(Q m ,t)/S(Q m ) clearly in- 
dicates that this length scale corresponds to the density 
fluctuation that gives rise to the main peak of S(Q). This 
length scale, 27r/Q m , is the size of the first shell of neigh- 
bours, the dominant clement of the liquid structure. The 
maximum of the relaxation time of F(Q, t) at Q m (De 
Gennes narrowing [15| ) is caused by the so-called "cage 
effect" U |23| whereby the structural relaxation is facili- 
tated by the diffusive motions of the particles within the 
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confining cages of the nearest neighbours. The cage size 
thus controls the relation between the diffusion and the 
structural relaxation, and its (^-invariance explains the 
latter's universality. 

The relaxation time for the Q-component of F(Q,t) 
scales as cx S{Q)/Q 2 0, 0> Ell Its variation with 
respect to Q and ip is shown in Fig. 21 It features a 
pronounced maximum at Q m for ip — 0.35 where the 
relaxation and diffusion follow the universal liquid rela- 
tion, see Fig. For tp = 0.2, that maximum disappears, 
and the decay of F(Q, t) is dominated by the density 
fluctuations with length scales exceeding 2n/Q m . The 
respective increase in the range of the diffusive motions 
dissipating these fluctuations results in the retardation of 
the relaxation process relative to diffusion as compared 
with the universal liquid behaviour, see Fig. [3] More- 
over, the stretching of &d(t) beyond a single-exponential 
behaviour can be interpreted as a superposition of mul- 
tiple exponential relaxation processes with a wide range 
of relaxation rates. 

We thus conclude that the observed universality of the 
liquid relaxation behaviour relative to diffusion is ex- 
plained by the confinement of the relevant density fluctu- 
ations to the first coordination shell which manifests itself 
as the maximum of the Q-dependent relaxation time at 
Q m (de Gennes narowing). 

We support this picture of liquid dynamics by the fol- 
lowing model. Consider a tessellation of space into unit- 
volume cubic cells. Each cell contains one point particle 
performing Fickian diffusion with the self-diffusion coef- 
ficient D. The diffusion starts from the centre of the re- 
spective cube, and it is confined to the latter by reflective 
boundaries. This dynamics relaxes the initial lattice con- 
figuration to the equilibrium perfect gas where the den- 
sity fluctuations are constrained by the single-occupancy 
confinement. This is supposed to model the confinement 
of the density fluctuations to the first neighbour shell 
that was concluded for the liquid dynamics. The relax- 
ation of the initial structure is assessed in terms of d(t). 
We note that, in this model, reflecting boxes' boundaries 
are equivalent to applying the d-minimizing permutation, 
and, for the uniform equilibrium distribution of particles' 
density within the cube, di = 1/4. The solution for d(t) 
can be obtained as the following serial expansion: 
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The decay of with respect to (r 2 (t)) = 6Dt shown 

in Fig. |21 demonstrates good agreement with the liquid 
relaxation data. 

These results highlight two principle aspects of the 
liquid dynamics representing essential conditions for the 
universal relation between the relaxation and diffusion, 
and the Stokes-Einstein relation. First, the relaxation 
process is dominated by the dissipation of the density 



fluctuations confined to the first coordination shell. This 
behaviour, also observed as de Gennes narrowing, can 
be diagnosed from the structure data as a maximum at 
Q m in S(Q)/Q 2 . Second, the trajectory of a diffusing 
particle facilitating the relaxation process is expected to 
be fully randomized within the first-neighbour distance, 
which implies Gaussianity of the diffusion process. The 
deviation from the universal pattern of liquid dynamics 
observed in the low-density domain is evidently caused by 
the presence of long-range density fluctuations violating 
the first of the above conditions. 
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